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'sj" I We give an analytic formula for the double distribution of hadronic in- 
variant mass and charged lepton energy to one-loop order of the perturbative 

, QCD. Although infrared singular, this quantity is closely related to physical 

^\ • observables that can be obtained thereof through proper convolution. 

^ : PACS numbers: 12.38.Bx, 13.20.He 

Oh. 1. Introduction 

(D 

1^ ' Today, the interest in bottom meson decays has never been stronger. Of 

those, the semileptonic kind are particularly fit for precision measurements 
of the Standard Model parameters such as the Cabibbo-Kobayashi-Maskawa 
matrix elements \Vcb\ and \Vub\ [1; 2] or the quark masses. The theoretical 
5^ \ predictions rely on perturbation theory as well as the effective theory of 

heavy quarks [3, 4, 5, 6] for inclusion of QCD effects. The present tech- 
niques are able to provide successful description of the semileptonic decays 
and it has been known that significant contributions to the B decays come 
from both perturbative [7, 8, 9, 10, 11] and nonperturbative [12, 13, 14] 
corrections. 

In this paper, we give the one- loop perturbative correction to the double 
differential distribution in terms of the hadron system mass and the charged 
lepton energy. The masses of the final particles are set to zero throughout 
the present calculation. 

The double differential distribution contains an explicit dependence on 
the fictitious gluon mass as it must due to the infrared singularities present 
both in the virtual and the real radiation corrections. It is not until after 
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one integrates over the invariant mass of hadrons that the apparent diver- 
gence cancels to leave a finite, measurable result. The distribution we give 
can be integrated over any desired range of the hadronic mass to yield an 
experimentally valid quantity. Hence, arbitrary cuts on the hadronic mass 
may easily be imposed. 

The result we obtain is consistent with the published formulae for the 
single differential distribution in both variables considered. In particular, we 
have checked they agree with the charged lepton energy distribution, first 
published by Jezabek and Kiihn[7], and the hadronic mass distribution, as 
found in [11, 15]. 



2. Kinematics 

2.1. Variables 

In our treatment of the decay of a i? meson, we need to consider both 
three- and four-body final states. The latter includes a gluon radiated by 
the initial or final quark. Denote the four-momenta of the involved particles 
as Q,q,l,i',G for the h quark, final state quark, electron, neutrino and 
gluon, respectively. Also define the hadronic four-momentum P = q + G. 
We are working in the rest frame of the B meson, which in the present 
approximation reduces to that of the b quark. 

The variables we use are x for the electron energy, t denoting the invari- 
ant mass of the leptons, and z for the invariant mass of the hadron system, 
all scaled according to the following formulae: 

X 2' 2' 2' V / 

with mf, the bottom quark mass. Also we will use the scaled masses of the 
gluon and the final u quark, which both pop up in the intermediate formulae 
in the course of calculating the contributions from three- and four-body 
processes. 

ma ^ ruu 

Xg = , e = ^/p= . (2) 

rub rub 

The system of the quark in the final state and the real gluon is described 

by the following quantities: 

P, = ^{l-t + z) (3) 



Ps = ^JPi-z = ^[l+t^ + z^-2{t + z + tz)]'^, (4) 
P±{z) = Po{z)±P3{z), (5) 
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where Po{z) and P'i{z) are the energy and the length of the momentum 
vector of the system in the h quark rest frame, and yp{z) is the corresponding 
rapidity. Similarly for the virtual boson W: 

Wo{z) = \{l + t-z), (7) 

W^{z) = ^W^ - t = ^[l + t^ + z^-2{t + z + tz)]^/^, (8) 

W±{z) = Wo{z) ± Wsiz), (9) 

Kinematically, the three body decay is a special case of the four body one, 
with the four-momentum of the gluon set to zero, thus resulting in simply 
replacing z = p. The following variables are then useful: 

PO = Po{p) = ^(1 - i + p) , P3 = P3{P) = ^JpI - P, (11) 
P± = P±{P) = PO ± P3 , W± = W±{p) = l- Pqr, (12) 

Yp = yj,{p) = \\n^-± , Y^ = yM = h^—- (13) 

^ P — ^ IX) — 

2.2. Kinematical boundaries 
The electron energy may vary within the range of 

0<x<l, (14) 

while for a fixed value of x, the allowed interval of z is 

0<z<l-x. (15) 

On the other hand, without fixing x, the hadronic mass may vary as 

0<2;<1, (16) 

and then the following boundaries on x are found for a specified z: 

0<x<l-z. (17) 

Given values of x and z, the intermediate boson invariant mass t varies 
within the following limits, 

0<t<x{l--^)=tm. (18) 
1 — X 
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3. Evaluation of the QCD corrections 

The differential decay rate for the process h ^ ulvi can be written as 

dV = dVo + dri,3 + dri,4 , (19) 

where 

dFo = Glml\VubfMi,3dn3{Q; q, I, i^)/7r^ (20) 
stands for the Born approximation, 

dri,3 = lasGWb\Vub\^Mi,3dn3{Q; q, I, (21) 

comes from the interference between the virtual gluon and Born amplitudes, 
and 

dri,4 = '^asG%ml\VubfMi,4dn4{Q; q, G, I, u)^ (22) 

is due to the real gluon emission. The Lorentz invariant n-body phase space 
is defined as 

dnn{P;pi,...,Pn) = 6^^\P -J2p^)1[.^. (23) 

3.1. Three-body contributions 

The Born approximation is evaluated according to (20) with the matrix 
element 

Mo,3 = q-lQ-i^, (24) 
and the Dalitz parametrization of the phase space, 

dH3{Q; q,l,u) = —dx dt. (25) 

Subsequently, the dot products appearing in Mo,3 are put in terms of the 
variables x and t, 

Q.q = + l.q = ^^{x - 1) 

Q-u = ^{1-p-x + t), l-v = ^t (26) 

Q-l = -^x, v-q = -^{1 — p — x). 

which readily yields the double differential distribution. On integration over 
t the latter yields the Born approximation to the electron energy spectrum 
(see Eqs. (39,41)). 
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The virtual correction is calculated along the same lines, now with the 
matrix element 



A1i,3 = - H^q ■IQ-U + H+pQ ■uQ-l + H^q ■ uq ■ I 
+^piH+ + H^)u ■ I + ^p{H+ -H_+ Hl)[1 -{Q-q-vm- v) 



1 



HL[l-{Q-q-v)]{q-v) , (27) 



substituted into (21). In (27), the form factors are as follows, 



Ho = 4(1 - Vo/Pa) In Ag + (2po/P3)[Li2 (1 - 



P-W- , 
P+W+' 



W- 



-Li2 (1 - ^ - YpiYp + 1) + 2(ln VP + YpWn, + Yp)] 

+[2p3Yp + {l-p-2t) In ^/t + 4, (28) 

H± = ^[l±il- p)/t]Yp/ps ± ^ In v^, (29) 

Hl = j{l - In Vp) + ^In Vp+ + f ^- (30) 



Integration over the invariant mass of the intermediate W boson gives the 
desired contribution to the double differential distribution, 



1 



34 



12ro dxdz 



=6{z){-log Xg[10x - 25x^ + —x-^ + (10 - 24x + 18x^ - Ax"^) log 



X 



-6x^ log e + 4.T''^ log e] + -(2x^ - 3x^)U2 (x) + — logx[121 

3 18 

-276x + 195x2 - 40x3 - log,T(30 - 72x + 54x2 - 12x3)] 

121 443 , 128 o 1 , „ , o 
+—X - —x^ + — x^ + - log e(-3x2 + 2x3 + Gx^ log e 
18 36 27 6 



-4x3 log e)}, 
where we have denoted 



X = 1 — X. 



(31) 



(32) 



While we present the decay rate assuming masslcss final particles, it is not 
possible to eliminate the final quark as well as gluon mass dependence out of 
the matrix element alone. This is removed once the integrated real radiation 
term is added. 
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3.2. Four-body contributions 
The four-body phase-space is decomposed as follows: 

dn4{Q;q,G,l,iy) = dz dUsiQ; P,l,u) dn2{P;q,G). (33) 

Discarding terms vanishing along with the gluon mass Xq — ^ 0, the matrix 
element representing the real gluon radiation may be written as 

^ _ _Bi B2 Bj 

''^ {Q-Gf QGPG^ {P- Gf ' ^ ' 

where 

B^ = q.l[Q.v{Q.G-l) + G-u-Q-uQ-G], (35) 
B2=q-l[G-u-q-vQ-G + Q-v{q-G-Q-G-2q-Q)] 

+Q-iy{Q-lq-G-G-lq-Q), (36) 
B3 = Q-u{G-lq-G- pl-P). (37) 

We refer the reader to [9] for the details of the integration that leads to the 
triple differential rate in terms oi x,t and z. This rate is conveniently split 
into terms according as they are infrared convergent or divergent. Hence 
we can write, 

^ = J'conv + Fdiv ■ (38) 

UTodxdtdz ^ ^ 

This expression is subsequently integrated over t. The infrared finite part 
has been integrated with the help of FORM. The formulae obtained in this 
way suffer from infrared and coUinear divergences. Thence the gluon and 
final quark masses subsist as regulators in spite of the limit we have taken. 
Of course, both remnant dependences vanish after the integration over the 
hadronic system mass is performed, which involves summing the virtual and 
real contributions. 



4. Results 

The first order QCD corrected double differential decay rate can be 
written in the form, 

= fo{x)5iz) + ^imz) + /f (x, z)eiz - X%)], (39) 



12ro dx dz Sir 

with 



^° - 192^' 
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where the first term on the right hand side is the Born approximation, given 
by 

(41) 



/o(x) = ^x\'i-2x), 



while 



-4/o(x) log^ \g + 7^(120x - 291x2 + 130x^)log Ac 
18 

+^(-10 + 24x- ISx^ + 4x^)log:rlog ( | 



A2 



+-(83 - 196x + 145a;2 - 32a;^) logx 



+-3-(249x - 426a;2 + 155x^ 
18 



^/o(x) [tt^ + 3Li2 (x)] , (42) 



and 



-[ - 2/o(a;) logz + ^ (l20x - 291^2 + 130x^) + ^(10 - 24x 
+18a;2 - 4x-^) log x] + [(lOz + Zz^)/x - [z + 2z^)lx^ 



2 2/^3ii^„ ^ I i;^-i(_692 - 25^2) + lx-2(ii^ + 252;2) 
2 2 ' 4 ^ ' 



3 X 
+x-^(-29zVl8) + (101/2 - VIxz - 81x + 45xV2 
-31z - 3zV2) log X + (32 - 6xz - 30x + ^x^ + 29z 



z 



16x^ 



+z ) log x log — + log ^(21x2: + 32x 

X 

-20xz + lOlx/2 - 97xV4 + 127z/4 + 283^736. 



9z - 5^2/3) 



(43) 



The above formula is easily integrated over either of the variables to give the 
single differential distributions in hadronic system mass or charged lepton 
energy. Then expressions confirming previous calculations [7, 11, 15] are 
found. The evident singularity of this distribution disappears after integra- 
tion over the hadronic system mass. That this indeed is so, can be seen by 
expressing it in terms of the following distributions. 



log z 



\ogz 



lim 
A->-0 



e{z-x) + -iog'xd{z) 



(44) 
(45) 



The substitution of these to Eqs. (42,43) results in the formal identification, 

Q{z -\l)\ = e{z - \l) - 5{z) log \l, (46) 
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^) -\W>^l5{z). (47) 



Upon application of Eqs. (46,47) to the correction terms, the latter take on 
the following form: 



1 



(-10 + 24j; - 18a;2 + 4a;^) log^ x 



+-(83 - 196x + Ubx^ - Z2x^) \ogx 

+-3-(249x-426a;2 + 155x3) -|/o(x)[7r2 + 3Li2(x) , (48) 

lo o L 



/f(x,z) = -2/o(x)(i^)^ + Q)^[^(l20x- 291x2 + 130x=^) 

+(^ - 8x + Gx^ - ^x^) log x] + [(lOz + 2,z^)/x - {z + 2z^)/x'^ 

+\z^/x^] log ^ + ^x-^(-69z - 25^2) + \x-'^{llz + 25^^) 

+x-3(-29zVl8) + (101/2 - 17xz - 81x + 45xV2 
-31z - Zz^/2) log X + (32 - 6xz - 30x + Ax'^ + 29^; 

+z^) log X log I + log z{2lxz + 32x - IGx^ - 9z - 52^/3) 

-20xz + lOlx/2 - 97xV4 + 127z/4 + 2832^36. (49) 

Clearly, the gluon mass does not enter the integrated distribution, de- 
fined as 

for which we obtain, 

Fi(x, z) = (Ci+C2 Z+C^ Z^+Ci Z^) logZ+C5 log^ Z+Ce-Z+CyZ^+Cg-Z^+Cg. (51) 

The coefficients ci to cg are as follows. 



ci = — [(120 - 288x + 216x2 - 48x3) log x + 120x - 291x2 + 130x3] , 

C2 = (-30x + 4x2 + 32)logx + 32x-16x2, 

C3 = 5x-^-^x"2 + ^(-6x + 29)logx+(21x-9)/2, 

C4 = x"^ - ^x~2 + ^x~3 + 1 log X - 5/9, 
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C5 = -x^/2 + x^/'i, 

eg = ^(-102x + 37x2 ^ - _^ ^3Q^ _ 4^2 _ 32) iog2 ^ _^ 

(74x - 33x2 )/4, 

C7 = ^(-28x - 91)logx + ^(6x - 29)log2x- 10x~Mogx- ^x~^ + 

—2/, - 13, 61 145 
X (logx + — ) - — x + — , 

cg = - 22 log X - 12 log2 X - 72x~^ log x - 162x~^ 

+48x~2 log X + 83x~2 - 16x~^ log x - 22x"^ + 101}, 

eg = —{6(83 - 196x + 145x2 - 32x^) log x + 12(-10 + 24x 
36 

-18x2 + 4x^) log2 X - 127r2x2 + 87r2x^ + 498x 

-852x2 + 310x3 - 36x2Li2 (x) + 24x3Li2 (x)}. (52) 

One way of making comparison between the parton model predictions and 
the resonance ridden experimental data is to consider moments of distribu- 
tion. We define those as 

While the zeroth moment corresponds to the electron energy distribution 
itself, the singular part of the double distribution leaves no trace on the 
higher moments. In fact, they are then vanishing in the Born approximation. 
The first five moments are expressed in terms of the following functions, 

2a 

Mn{x) = -^m„(x), n > 1, (54) 

which are given by the formulae beneath (see also Fig. 1): 

mi = (-35/144 + 5x2/8 - 4x^9 + x^/16) log X 

-35x/144 + 19x2/72 -xV48, (55) 

m2 = (-449/3600 + 7x/24 - 61x2/360 - 2.X-V45 + 13x^/240 
-13x^1800) log X - 449x/3600 + 533x2/1800 - 9x^40 
+109x^/1800 - xVl44, (56) 

mg = (-103/1800 + 119x/600 - 29x2/120 + 19x^/180 + x^/120 
-3xV200 + xVeOO) log X - 103x/1800 + 697x2/3600 
-653x^2700 + 49x''/360 - 7xV200 + 47xVl0800, (57) 
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Fig. 1. Moments of hadron mass distribution function of the scaled electron 

energy x. The first five moments are shown as indicated by the legend. 

m4 = logx(-1313/44100 + 2x/15-81arV350 + 59x^315 -5x^/84 

-x^/175 + 2x^315 - 2x'^/3675) - 1313x/44100 + 5729xV44100 
-11x^/49 + 857x''/4410 - 781xV8820 + 311x^/14700 
-19x^/7350, (58) 

7715 = logx(-485/28224 + 23x/245 - 211x^1008 + 151x^630 

-95x^/672 + 2x^63 + 29x^/5040 - .t'^/294 + 31xVl41120) 
-485x/28224 + 1321.tV14400 - 15879x'V78400 
+33599xVl41120 - 4523xV28224 + 979x^15680 
-9809x^/705600 + 163xVl00800. (59) 



5. Summary 

One-loop QCD correction has been found to the electron energy and 
hadronic system mass double distribution in semileptonic B decays. The 
result has not been published before. It agrees with the known single dis- 
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tributions in both variables. Although infrared singular, it is of use for 
experimental analysis with appropriate cuts imposed. 
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